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We present Fermi-LAT and multi-frequency, multi-epoch VLBA data for the TeV blazar Mrk 421. We collected 
the data during a long and intensive multi-frequency campaign in 2011. We study the gamma-ray light curve, 
the photon index evolution and their connection to the radio data on sub-parsec scales, including total intensity, 
polarized flux density, polarization angle, spectral index, and rotation measure both in the core and the jet 
region. The VLBA data were obtained at 15 and 24 GHz for 12 epochs and at 43 GHz for 23 epochs, thus 
providing the best temporal and spatial coverage in the radio band ever achieved for a TeV blazar. We provide 
significant constraints on the jet Doppler factor, the presence of proper motion, the magnetic field configuration, 
and an intriguing connection between variability in the radio data and the gamma-ray light curve: the total 
intensity and polarized core emission reach a peak simultaneously to the main gamma-ray peak, followed by 
a rotation of the polarization angle at low frequency. Opacity-related, long wavelength polarization swings 
are also detected later in the year, possibly related to secondary peaks in the gamma-ray light curve, setting 
constraints on the physics of the gamma-ray zone. 


1. INTRODUCTION 

Multi-wavelength variability studies provide an ex¬ 
traordinary opportunity to break degeneracies be¬ 
tween the various blazar emission models [suniiin], 
which predict flux variations (at a given energy band) 
with particles of different energies, cooling times, and 
cross sections for different processes [7]. However, 
the complexity in resolving the underlying processes 
occurring in blazars can only be achieved through 
a well-sampled dedicated monitoring from radio to 
gamma rays lasting several years. Indeed, some of the 
latest and most interesting results on blazars come 
precisely from variability studies from well-sampled 
coordinated years-long multi-instrument observations 
[Diais]. Unluckily, these multi-instrument observa¬ 
tions were performed on flat spectrum radio quasars, 
but not on TeV BL Lac objects [with the exception 
of BL Lacertae, la, perhaps owing to the fact that 
these sources were challenging to observe in radio 
and gamma rays. The advent of Fermi-LAT and the 
current generation of Cherenkov telescopes (H.E.S.S., 
MAGIC and VERITAS) permits the accurate and sys¬ 
tematic (regardless of the activity level) monitoring of 
the high energy bump, where a large fraction of the 
blazar emission is produced. In addition, recent works 
have shown that TeV blazars can be successfully im¬ 
aged with the VLBA at up to 43 GHz, revealing the 
core and the inner jet structure with great accuracy 
[14j . Therefore, the currently available instrumenta¬ 
tion allows to study these sources with an unprece¬ 
dented level of detail, which has the potential to shed 
light on the understanding of these complex objects. 

In order to address the challenge involved in break¬ 
ing the degeneracy between models, we have orga¬ 


nized the most ambitious multifrequency monitoring 
to date of the famous TeV BL Lac object Mrk 421, 
covering sub-mm (SMA), optical/IR (GASP), UV/X- 
ray {Swift, RXTE, MAXI), and gamma rays {Fermi- 
LAT, MAGIC, VERITAS). We summarize here the 
main results of the multi-frequency, multi-epoch, full 
polarization VLBA observations and of the gamma- 
ray Fermi-LAT monitoring through 2011 [aiiiio]. 

Throughout this paper, we adopt Hq = 70 km sec“^ 
Mpc“^, LIm = 0.25 and Ha = 0.75, in a flat Uni¬ 
verse, so that at the redshift of Mrk 421 {z = 0.031), 
Imas = 0.59 pc. Spectral and photon indexes a and 
r are defined such that the radio flux density and the 
gamma-ray photon flux are proportional to and 
E~^, respectively; angles are measured from north 
through east. 


2. OBSERVATIONS 
2.1. Radio observations 

We observed Mrk 421 with the VLBA for 12 times 
(once per month) throughout 2011 at 15, 24, and 
43 GHz; at the latter frequency, we expanded our 
dataset with 11 more observations obtained through 
the VLBA-BU-BLAZAR monitoring program. We 
carried out a full calibration and analysis describing 
the evolution with time of total intensity and polar¬ 
ized flux density, and of their combinations such as 
spectral index and rotation measure. Core and in¬ 
ner jet are spatially resolved and separately analyzed. 
Spatial resolution is as good as 0.2 pc at 43 GHz and 
sensitivity just below 1 mJy/beam. 
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Figure 2: Mrk421 43 GHz core parameters as a function 
of time. From top to bottom panel: gamma-ray light 
curve (0.1 to 100 GeV, in 10“® ph cm“^ s“^ units, shown 
for comparison), total intensity flnx density, polarized 
flnx density, fractional polarization, EVPA. In the total 
flnx density panel (second panel from top), filled circles 
indicate data from the main set of observations, empty 
triangles indicate additional data from the 
VLBA-BU-BLAZAR monitoring project. 


2.2. Gamma-ray observations 

We analyzed gamma-ray data from the Large Area 
Telescope on board Fermi, which continuously scans 
the whole sky in the energy range 100 MeV < E <300 
GeV. We analyzed the data with the ScienceTools 
software package version v9r32p5, using instrument 
response functions P7REP_SOURCE_V15 and follow¬ 
ing the standard procedures^. Mrk421 is bright 
enough to be significantly detected in every weekly 
bin, and we obtained photon flux and photon index 
for every bin. 


3. RESULTS 


We present sample images at the three frequencies 
in total intensity and polarization in Figure In 
the following subsections, we describe the main results 
about the radio emission from the core (Sect. [3^ and 
jet (Sect. 3.2), and the unresolved gamma-ray source 


(Sect. 3.3). 


3.1. Radio core 

The core is the most prominent feature in the radio 
images. It is bright and compact at all epochs and 


Figure 1: VLBA images of Mrk421 at 15 GHz (top, 2011 ^httpiZ/ferml.gsfc.nasa.gov/ssc/data/analysis/ 

January), 24 GHz (middle, 2011 February), and 43 GHz documentation/Cicerone/ 

(bottom, 2011 April). 
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Figure 3: Mrk421 core rotation measure (top panel) and 
intrinsic EVPA (bottom panel) as a function of time 
throughout 2011. 


JET 15GHz 



Figure 4: Mrk421 15 GHz jet parameters as a function 
of time. From top to bottom panel: total intensity flux 
density, polarized flux density, fractional polarization, 
EVPA. 


frequencies and it shows variability in total intensity 
and polarization. We show its light curve and the time 
evolution of its polarization properties in Figure 
along with the gamma-ray light curve that is described 
later on (Sect. 3.3). 


The radio core light curve showed a broad peak 
around 2011 February-March at all frequencies. This 
peak is particularly prominent at 43 GHz, thanks to 
the improved sampling offered by the additional 11 
observations from the Boston University blazar mon¬ 
itoring project. In particular, the core flux density 
reached its peak of S '43 = 415 mJy on 2011 March 1 
at 43 GHz. The core spectrum is generally flat, with 
average spectral indexes = 0.16 and = 0.43. 
The spectral index shows a flatter when brighter be¬ 
haviour. 


The core polarization properties are also variable: 
the polarization fraction varies as a function of both 
frequency (being ~ 1% at 15 and 24 GHz and ~ 2% at 
43 GHz) and time, with a (3.6 ± 0.5)% peak on 2011 
March 29. The electric vector position angle (EVPA) 
X has some fluctuations, which are more prominent at 
low than at higher frequency. At 43 GHz, it oscillates 
slightly with a mean value of 144° and a standard 
deviation of 17°, being overall well aligned to the jet 
axis. At 15 GHz, we detect two prominent 90° flips in 
July and September, most likely due to opacity; even 
after correcting these values, there is larger variability 
with a mean of 119° and a standard deviation of 29°. 

Epoch by epoch, we performed a x vs fit to de¬ 
termine the rotation measure RM and the intrinsic 
EVPA xo in the source, under the assumption that 
the frequency dependence of the observed EVPA’s is 
due to Faraday rotation effect. The evolution of RM 
and Xo is shown in Figure The RM values oscillate 
between —3000 and 2000 rad m“^. 


3.2. Radio jet 

The VLBA data clearly detect a one-sided jet at all 
three frequencies. The jet extends for about 5 mas 
(3 pc) in PA —35° at 15 GHz and somewhat less as 
frequency increases. In Figure we show the total 
intensity and polarized flux density, the fractional po¬ 
larization, and the EVPA of the jet emission at 15 
GHz, from top to bottom. 

We do not detect any prominent variation in neither 
the jet’s structure nor its total intensity flux density. 
In particular, we did not detect any significant super¬ 
luminal motion of components laiin]; as a matter of 
fact, we do not detect well defined, compact jet com¬ 
ponents at all. At a reference distance of 1.5 mas from 
the core, the jet-counterjet brightness ratio is i? > 30. 
The spectrum is steeper than in the core. 

In polarization, we detect emission from the jet at 
15 and 24 GHz, with a fractional polarization value 
that is much larger than in the core (^ 15%). How¬ 
ever, given the steep spectrum of the jet and its low 
surface brightness, at 43 GHz the polarized emission is 
below the noise level. The EVPA in the jet is remark¬ 
ably stable, with a mean of 61° and a standard devi¬ 
ation of just 9°. Given the jet axis PA of —35°, this 
means that the EVPA is nearly orthogonal to the di¬ 
rection of the flow, and the magnetic field well aligned 
with it. 


3.3. Gamma-ray source 

Mrk421 is point like at the Fermi angular reso¬ 
lution. This is not only true for the angular scales 
probed by the VLBA observations presented here but 
also when we consider the total extent of the source 
as mapped by compact interferometers best suited to 
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Figure 5: Mrk421 gamma-ray parameters as a function 
of time. Top: photon flux; bottom: photon index. The 
dashed lines indicate the 1-yr mean value. 

map a possible large scale diffuse emission [^. How¬ 
ever, the variability time scales observed for the high 
energy emission and the broadband spectral energy 
distribution modelling [5] clearly suggest that most of 
the gamma-ray photons are produced in a compact 
region of size smaller than the radio core itself. 

During 2011, the mean photon flux and the pho¬ 
ton index are F = (17.4 ± 0.5) x 10“® ph cm“^ s“^ 
and r = 1.77 ± 0.02, respectively. The light curve 
on weekly time bins shown in Figure (top panel) 
reveals variability, with three peaks in 2011 March, 
September, and November. The brightest week is 
that between 2011 March 5 and 11, (MJD 55 625- 
55 631), with a gamma-ray flux as large as (38± 11) x 
10“® ph cm“^ s“^. The photon index on the other 
hand is quite stable and we do not find any evidence 
for spectral variability (see Figure]^ bottom panel). 


4. DISCUSSION 

The results from the radio observations suggest that 
the jet is not strongly Doppler boosted already on par¬ 
sec scale. This stems from the low brightness temper¬ 
ature estimated both from the variability time scales 
(T'B,var 2.1 X 10^° K) and the core size and flux 
density (Tb ~ 10^^ K) [TU]. On the contrary, TeV 
variability on time scales of about 30 minutes [1] , and 
both hadronic (synchrotron proton blazar) or leptonic 
(one-zone synchrotron self-Compton scenario) fits to 
the broadband spectral energy distribution require 
Doppler factors 6 in excess of ^ 10, and possibly as 



Figure 6: DCF results for Mrk421, considering 
gamma-ray Fermi data and 43 GHz VLBA data. Black 
line: observed data; gray curves: 99.7% confidence level 
threshold obtained from the combination of different 
PSD slopes. 

large as <5 = 50 |2]. This so-called “Doppler factor 
crisis” characterizes also other TeV blazars; the most 
natural way to solve this crisis is to localize the radio 
and the gamma-ray emission zones in different regions; 
typically, a velocity structure of the jet is assumed, ei¬ 
ther along or across the jet axis. 

In our case, we have the possibility to compare the 
radio and gamma-ray light curves and to carry out a 
discrete cross-correlation function (DCF) analysis for 
the two datasets. In particular, we carried out a DCF 
analysis over a range of radio-gamma ray delays be¬ 
tween — 100 and -1-100 days, with a bin of 15 days. We 
find a peak for the correlation (0.54) is obtained for 
At = 0 day delay. In order to assess the significance 
of this peak, we generated 3000 fake light curves with 
the same mean and standard deviation as the observed 
light curves but with variable power spectral densities 
(PSD). We show in Figure]^ the results of our DCF 
analysis on the real data compared to those obtained 
for various combinations of the simulated ones. In 
particular, the peak is significant at the > 99.7% con¬ 
fidence level for the combinations of light curves with 
short time scale gamma-ray variability (/3-y_ray < 1-5) 
and longer time scale radio variability (1 < /3r < 2.5). 

As for the magnetic field configuration in the radio 
core and jet, we found somewhat discrepant results. 
The core B-field perpendicular to the jet axis is typ¬ 
ical of shocked regions; in the jet, the magnetic field 
is aligned with the jet axis; possible interpretations 
of this result include stretching of an initially trans¬ 
verse magnetic field by a layered velocity structure 
and an helical field with a small pitch angle. Finally, 
we note that the variability features in the polarized 
flux light curve indicate also a connection between 
magnetic field and gamma-ray emission. 
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